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ABSTRACT 

We explore a new possibility of electroweak baryogenesis in the next-to-minimal supersym- 
metric standard model. In this model, a strong first-order electroweak phase transition can 
be achieved due to the additional singlet Higgs field. The new impact of its superpart- 
ner (singlino) on the baryon asymmetry is investigated by employing the closed-time-path 
formalism. We find that the CP violating source term fueled by the singlino could be 
large enough to generate the observed baryon asymmetry of the Universe without any con- 
flicts with the current constraints from the non-observation of the Thallium, neutron and 
Mercury electric dipole moments. 
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The origin of the baryon asymmetry of the Universe (BAU) is one of the longstand- 
ing problems in particle physics and cosmology The BAU is observed from the primordial 
abundances of the light elements (D, 3 He, 4 He, and 7 Li) induced by the big-bang nucle- 
osynthesis (BBN) and the cosmic microwave background (CMB) data. From these data 
the baryon-to-entropy ratio is found to be pQ 



The baryon asymmetry of the Universe must arise after the inflation if it happened, and 
before the BBN era. As pointed out by Sakharov j2] the following three conditions must 
be satisfied in order to generate the BAU: (i) baryon-number violation, (ii) C and CP 
violation, and (iii) the departure from thermal equilibrium. The last condition could be 
evaded if the CPT symmetry is violated. Although the standard model (SM) of particle 
physics can in principle satisfy all the above three conditions, it turns out that the CP 
violation coming from the Cabbibo-Kobayashi-Maskawa matrix [3] is too feeble to generate 
sufficient BAU [I]. Furthermore, the electroweak phase transition (EWPT) in the SM is a 
smooth crossover for a Higgs boson with a mass above the LEP bound (~ 115 GeV) [5], 
and therefore renders the condition of departure from the thermal equilibrium infeasible. 
Thus, physics beyond the SM is indispensable to address this issue. 

Supersymmetry (SUSY) is an attractive framework for new physics since it can pro- 
vide elegant explanations for a number of key questions that cannot be accommodated 
in the SM [6]. Supersymmetry can stabilize the electroweak-symmetry breaking scale 
by providing a natural cancellation mechanism of the quadratic divergences thus solv- 
ing the gauge-hierarchy problem, offer a good candidate for the cold dark matter, enable 
the gauge-coupling unification, break electroweak symmetry dynamically, and generate 
sufficient BAU. The minimal realization of SUSY, known as the minimal supersymmetric 
standard model (MSSM), is very attractive because of its simplicity. It is, however, known 
to suffer from the so-called /i problem [7] and the little-hierarchy problem as noted more 
recently. Moreover, the electroweak baryogenesis (EWBG) scenario [5] in the MSSM is 
highly restricted by current experimental data [9TTTT]. Especially, a Higgs boson which has 
escaped the detection at LEP 2 [12] and/or perhaps weighs around 125 GeV as hinted by 
the recent ATLAS/CMS results [13] has pushed the MSSM to the edge of the allowed pa- 
rameter space that can be consistent with a first-order EWPT [TTJ . This tension, however, 
could be relaxed if the MSSM is extended [HH2D]. 

In this Letter, we consider the next-to-minimal supersymmetric standard model 
(NMSSM) [21], in which the fi problem is solved naturally. It has been shown that a strong 
first-order EWPT can happen in the NMSSM much more easily than in the MSSM [T^|[T5]. 
In the presence of the additional gauge-singlet field, the NMSSM Higgs potential contains 
explicit/spontaneous CP violation even at the tree level [22] • Furthermore, this new ad- 
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ditional source of CP violation may give rise to sufficient BAU even when the MSSM 
fails. 

The NMSSM superpotential is 
H^nmssm = U c h u QH u + D c h d H d Q + E°h e H d L + XSH u H d + | S 3 , (2) 

where S denotes the singlet Higgs superfield, H ud are the two SU(2) L doublet Higgs su- 
perfields, and Q, L and U c , D c , E c are the matter doublet and singlet superfields, re- 
spectively, related to the up- and down-type quarks and charged leptons. We note that the 
superpotential respects an extra discrete Z3 symmetry. 

The tree-level Higgs potential includes the terms coming from the soft SUSY-breaking 
terms: 

Koft = m\H\H d + m 2 2 HlH u + m 2 s \S\ 2 + (\A x SH u H d - ^kA k S 3 + h.c.J , (3) 

in which A, k, A x , and A K may contain nontrivial CP phases. After the neutral components 
of the two Higgs doublets and the singlet develop their vacuum-expectation- values (VEVs), 
one may have three rephasing-invariant combinations of the CP-odd phases in the tree-level 
potential [22] : 

& + 0a*; <P'x + <I>a k (4) 

where <p' x = <p\ + 9 + tp and <p' K = cf) K + 3(p with 9 and (p parameterizing the overall CP 
phases of the doublet H u and the singlet S, respectively. It turns out that the latter two 
CP phases are fixed up to a two-fold ambiguity by the two CP-odd tadpole conditions, 
while the difference <ft' x — <p' K remains as an independent physical CP phase. This tree-level 
CP violation is a distinctive feature of the NMSSM compared to the MSSM. 

The stringent constraint on the tree-level CP phase <p' x — <p' K may come from the 
non-observation of the electric dipole moments (EDMs) for Thallium [23] , the neutron [23] , 
and Mercury ( 199 Hg) [251 EE]. The detailed study shows that the maximal CP phase 
<p' x — cf)' K ~ 90° could still be compatible with the current EDM constraints taking account 
of the uncertainties in the calculations of the Mercury EDM when the sfermions of the first 
two generations are heavier than about 300 GeV [27] . 

In this Letter, we address the question whether sufficient BAU can be generated by 
the CP phase <p' x — <p' K ~ ±90° when the CP phases appearing in all the other soft SUSY- 
breaking terms are not large enough for the BAU because of the tight EDM constraints with 
the SUSY particles within the reach of the LHC. Explicitly, we are taking sin(<// A + <pA f ) = 
sin(0' A + <pi) = with <pA f and <pi denoting the CP phases of the soft trilinear parameters 
Af and the three gaugino mass parameters M i=li 2,3, respectively. 
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The scenario we are considering has an intermediate value of tan (3 with small vg ~ 
v(T = 0) ~ 246 GeV: 



tan/3 = 5, w 5 = 200 GeV, 

|A|=0.81, |«|=0.08; |A A | = 575 GeV , = 110 GeV ; 

<j>' x - <f>' K = ±90° , sign [cos(0' K + A J] = sign [cos(0' A + <f> Ax )] = +1 . (5) 

The other parameters are chosen as 

Mq 3 = 1 TeV , il% 3 = 150 GeV - 1 TeV , M Bs = 250 GeV - 1 TeV ; \A t \ = \A b \ = 1 TeV . (6) 

We find Ma — 600 GeV for the parameters chosen and have taken Mi = M 2 = —200 
GeV to fix the neutralino and chargino sectors. We find that a first-order phase transition 
could occur in some regions of the parameter space of this scenario which is needed for the 
EWBG. Actually, the diverse patterns of the EWPT in the NMSSM have been investigated 
in Ref. [15J. Among the patterns the so-called type-B transition opens a new possibility for 
the strong first-order EWPT. In this type of transitions, the lighter stop could be heavier 
than the top quark, in contrast to the MSSM EWPT. Instead, the singlet Higgs field plays 
an essential role. In a typical parameter-space point, during the EWPT the VEVs (v, Vg) 
change from (0,600 GeV) to (208 GeV, 249 GeV) at the critical temperature T c = 110 
GeV. The dramatic change of vg is the most important feature of the type-B transition. 
As discussed in Ref. [15], the relatively small k is required to ensure the local minimum in 
the vg direction as in the scenario Eq. (JSJ). 

For the calculation of the baryon density, the Closed-Time-Path (CTP) formalism is 
employed [2S],[29] and we closely follow Refs. [221120] B- The baryon density in the broken 
phase is given by 



Kb = — - — A 



r 1 + r 2 ^L U-SX 
T„D \ D 



v w [Dv w + (2D q - D)^Jvl + A7lD q \ + A7lDD q ' 

where rip = 3 is the number of fermion families, T ws = 6 k a\ T ~ 0.5 x 10 -5 T with k ~ 20 
and «2 — 1/30, and T ss = 16k' a^T with k' = 0(1). We are taking the bubble wall velocity 
v w = 0.04. The relaxation term 1Z = T ws [|(1 + nsg/6)^ 1 + |] with n sq for the number of 
light squark flavors. The diffusion constant is given by 

D = {9k Q k T + k B k Q + Ak T k B )D q + k H {9k T + 9k Q + k B )D h 
9k Q k T + k B k Q + Ak T k B + k H (9k T + 9k Q + k B ) 



*For the more precise calculation, the profiles of the bubble wall and the profile-dependent masses and 
widths should be taken into account when solving the quantum transport equations during the EWPT. 
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with D q = 6/T and Dh = 110/T ^T\. The k factors are given by the sum Ivq^t,b = 
k QL,t R ,b R + k^ L t R i R and k H = k Hd + k Hu + k & where 

6 f°° e x I 

k = 9i— dx . 2 xJ x 2 - m 2 /T 2 . (9) 

Note Qi = 1 for a chiral fermion and a complex scalar. For example, g^ti) = gi(ti) = 3 
taking account of the 3 colors, gi{H) = 2 for a Higgs doublet, and gi(H) = 2 for a Dirac 
Higgsino. We further note that ki = gi and ki = 2gi for a chiral fermion and a complex 
scalar, respectively, in the zero- mass limit. The coefficients r 12 are given by 

r = 9k Q k T - 5k Q k B - 8k T k B ^ = k 2 B (5k Q + 4k T ){k Q + 2k T ) 

n ~ k H {9k Q + 9k T + k B ) ' r2 ~ k H {9k Q + 9k T + k B ) 2 ' 1 j 

Note that, in the limit of very heavy stops and sbottoms, the sfermion contributions to 
their corresponding k factors become negligible and one may have /cq = 6 and kx = k B = 3, 
which leads to vanishing n to which the baryon density n B is directly proportional. We 
also observe that the coefficient r\ vanishes when the stops and sbottoms are all degenerate. 
In the calculation of the k factors, we have taken the thermal masses given in Ref. [32] with 
appropriate modifications. Finally, the parameter A is given by 




A^k H L w J^ , SH ° , (111 



gCPV 

i>HL w \ — — ; 

where L w denotes the bubble wall width and rp = F/(r^ + T/J. The CP-conserving 
particle number changing rate is 



{9k Q + 9k T + k B ){T M + Y h 



(12) 



" 9k Q k T + k B k Q + Ak T k B + k H (9k T + 9k Q + k B ) ' 
where 

6 / \ 6 

— 7p l^t + ) ' ^ h = (^H ± W ± + ^H°W° + ^5o B o + Ffios) , (13) 

with Tr}y = r~~ — r~~. We refer to Ref. |29] for details of the calculations of the rates 

HX HX HX 1 1 

F~~ and F~ 

t,t H±W± ,H°W° ,H°B° ,H°S~ 

Here we present the analytic expression for the singlino-driven CP-violating source 
term appearing in the parameter A in Eq. (Ilip as follows: 

= -2|A| 2 |M g ||/i cff |^/3sin(0 A - K )Z|~ O , (14) 
where |/x e g | = \\\v s /V2 and 



M § (T)\ 



2 1 V2 



|A| +2|K| T 2 



(15) 
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including the singlino thermal mass term. We assume that there is no spontaneous CP 
violation or 9 = if = 0. We note that the source term vanishes when sin(0A — <f> K ) = or 
$ — 0. If (3 has a kink-type profile, $ ~ v w Aj3 /L w and A becomes independent of L w . In 
the MSSM, Af3 was found to be O(10~ 2 - 1(T 3 ) We are taking A/3 = 0.02 in our 

estimation of the source term Q. The fermionic source function takes the generic form 
of 

k 2 



1, 



Air 2 



dk 



UJiUJj L 



(l - 2Re(nj))l(ui, F h uj, Vj) + (l - 2Re(nj)) /(wj, E,-, u iy r. 



2(lm(ni) + Im(nj)) G(wi, T i5 wy, r 3 -) 



where 



1 



11; 



(16) 



(17) 



with ojj = \fk 2 + mf. We note that the thermal width Tj at finite temperature is given 
by the imaginary part of its self-energy which is nonvanishing independently of whether 
the particle i is stable or not [34] • Specifically, for the calculation of the source function, 
we have taken the thermal widths given in Ref. [35] 0. The thermal functions I and G are 
defined by 



I(a, b, c,d) = - 



1 



1 

+ -x 



2[(a + c) 2 + (b + d) 2 ' 
1 



sin 



G(a, 6, c,d) = — 
1 

+ x 



2 [(a 
1 



(b 
1 



df 



sin 



2 arctan 



2 arctan 



a + c 



2 [(a + c) 2 + (6 + d) 2 ; 
1 



cos 



2 arctan 



b + d 
a + c 



2 [(a 



(b + 



cos 



2 arctan 



b + d 
a — c 



We note the thermal functions lead to oc 



Tj when ~ , 



6 + c/_ 

T>r 8 



(19) 



In Fig. (TJ we show our predictions for the singlino- driven Yb/Y^ as functions of 
v s {T c ) taking three values of M Bs = 250 GeV, 400 GeV, and 1 TeV with Y° h being the 
averaged value of the two BAU data in Eq. ([T]). Here vsiTc) denotes the singlino VEV at 
Tc and it can take on any value between 250 GeV and 600 GeV in the type-B EWPT: 



{v, v s ) = (0, 600 GeV) -> (208 GeV, 249 GeV) 



(20) 



' Note the source term grows linearly with A/3 and our choice is optimal. We observe that the variation 
depending on the choice of A/3 is to be regarded as the theoretical uncertainty and the precise determination 
of A/3 in the NMSSM is beyond the scope of this Letter. 

+For the thermal width of the singlino, we are taking = 0.03 T considering the large coupling 
|A| = 0.81. We find Yg is affected by the amount of about (25-35) % as Tg varies between 0.003 T and 
0.03 T. 
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Figure 1: The singlino- driven Yb as functions of the singlino VEV at the phase transition 
which, in principle, can take any value between 250 GeV and 600 GeV in the type-B EWPT. 
We are normalizing our predictions to Y^ h = 8.8 x 1CT 11 for three values of Mg 3 = 250 
GeV, 400 GeV, and 1 TeV We fix Mq s = Mq z = 1 TeV 



at Tq = 110 GeV [TBJ. We find that the singlino-Higgsino mass difference becomes smaller 
as vs decreases, leading to the larger Yg. The more accurate determination of Yg requires 
the knowledge of the profiles of the bubble wall and the treatment of the diffusion equation 
beyond the formalism developed in Refs. [291130]. We leave the more precise determination 
of Yb in the NMSSM framework for future work 



From Fig. [TJ we can see that Yb is much suppressed when Mg 3 = 1 TeV since the 
ratio ri is almost vanishing when both the stops and sbottoms are heavy and/or degenerate. 
However, Yb grows quickly as Mq 3 decreases. When M^ 3 = 400 GeV, the ratio Yb/Y^ 
is larger than 1 in the region vs{Tc) 440 GeV. Furthermore, in the case with M^ 3 = 
250 GeV, sufficient BAU can be generated via the singlino-driven mechanism, irrespective 
of the nonlinear dynamics during the EWPT. We found the similar behavior by fixing Mq 3 
and Mfj^ and varying M^ 3 , as shown in Fig. 2. We observe the same mass M^ 3 would 
give a slightly smaller Yb compared to the same mass of . Summarizing these results, 
we conclude that a sizable mass splitting either in the stop sector or the sbottom sector 
is needed for the successful singlino-driven BAU. In passing, we note that the resonance 



7 



100 



10 



M 5a = 1000 GeV 



J3 



C 1 



o.i 



o.oi 



Mr T = 150 GeV 




1 — I — | — I — I — I — I — | — I — I — I — I — | — |— 

200 300 400 500 

v s (T c ) [GeV] 



600 



700 



Figure 2: The same as in Fig. U\ but varying : Mfj 3 
and 1 TeV. We fix M 5 =M„ =1 TeV. 



150 GeV, 250 GeV, 400 GeV 



enhancement of the widths T^ ± ^ ± ^ ^o h°b° induces the dip around vsiTc) 
where \Mi\ = \M 2 \ = |// e ff|- 



350 GeV 



Lastly, the EWBG scenario considered here includes the two light Higgs states well 
below 100 GeV escaping LEP constraints [22] and the lightest neutralino of about 45 GeV 
with the singlino fraction of ~ 40 %. The light Higgs bosons and the lightest neutralino de- 
serve further studies in connection with the current LHC Higgs searches and the abundance 
of dark matter in the Universe 

In this Letter, we have examined a new possibility of a singlino-driven mechanism for 
the BAU in the NMSSM framework. In contrast to the MSSM, explicit and/or spontaneous 
CP violation can occur in the NMSSM Higgs potential even at the tree level. We emphasize 
that this new source of CP violation may solely give rise to the sufficient BAU without any 
conflicts with the current EDM constraints, as long as there is a sizable mass splitting in 
the stop and/or the sbottom sectors with the lighter stop and/or sbottom weighing below 
~ 500 GeV 
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